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We describe the relation between spin fluctuations and superconductivity in a highly-ordered 
sample of YBa2Cu306.5 using both polarized and unpolarized neutron inelastic scattering. The spin 
susceptibility in the superconducting phase exhibits one-dimensional incommensurate modulations 
at low-energies, consistent with hydrodynamic stripes. With increasing energy the susceptibility 
curves upward to a commensurate, intense, well-defined and asymmetric resonance at 33 meV with 
a precipitous high-energy cutoff. In the normal phase, which we show is gapless, the resonance 
remains surprisingly strong and persists clearly in Q scans and energy scans. Its similar asymmetric 
spectral form above Tc=59 K suggests that incoherent superconducting pairing fluctuations are 
present in the normal state. On cooling, the resonance and the stripe modulations grow in well 
above Tc below a temperature that is comparable to the pseudogap temperature where suppression 
occurs in local and low-momentum properties. The spectral weight that accrues to the resonance is 
largely acquired by transfer from suppressed low-energy fluctuations. We flnd the resonance to be 
isotropically polarized, consistent with a triplet carrying 2.6 % of the total spectral weight of the 
Cu spins in the planes. 

PACS numbers: 74.72.-h,75.25.+z,75.40.Gb 



I. INTRODUCTION 

Spin fluctuations play a fundamental role in the super- 
conductivity of high-temperature superconductors. In- 
tensive research has spawned a plethora of theories for 
the interplay between superconductivity and antiferro- 
magnetism but no model is generally accepted for the 
unusual behavior of the doped planar cuprate supercon- 
ductors. The search for new phases of matter, especially 
in the undcrdopcd region, as the boundary with the an- 
tiferromagnetic phase is approached, continues to pro- 
vide stimulus for experiment and theoryi Many recent 
theories for both the normal and superconducting states 
have made specific predictions for the spin correlations 
and therefore a detailed and complete study of x"(Qi i^) 
is essential to understanding the cuprates. 

One of the most interesting properties of the cuprate 
phase diagram is the existence of a pseudogap phe- 
nomenon seen clearly in tunneling, NMR, transport, and 
ARPES experiments!^ Tunneling data have evinced clear 
evidence of a gap-like structure in the density of states 
well above the onset to superconductivityi^ Some of the 
cleanest data have come from NMR studies which have 
found a suppression of both the Knight shift and the re- 
laxation rate 1/TiT in the normal statei^ 

Many theories have been constructed to explain the 
pseudogapi^ Most recently an orbital current or flux 
phase has been suggested which predicts circulating cur- 
rents flowing in the copper-oxide planesi2i2ii£ In the 
static version of the theorj*iS*ii a hidden order param- 
eter, called the d-density wave (DDW), condenses at a 
temperature well above the superconducting transition 



temperature. Since, the static DDW order would pro- 
duce a moment, this phase should be directly observable 
with neutrons as new Bragg peaks, and since orbital cur- 
rents break an Ising-like symmetry, a gap should exist in 
the spin excitation spectrum. We have previously found 
that no new Bragg peaks are observable in the normal 
phase of the ortho-II ordered YBa2Cu306.5 (YBCOg.s), 
thus excluding strong DDW order that breaks the sym- 
metry at Q = (7r,7r). We cannot exclude the possibility of 
broken orbital symmetry at zero wave vector— or the pos- 
sibility that the orbital currents are primarily dynamic^ 

Another theory to explain the variety of phases 
throughout the cuprate phase diagram is the stripe the- 
ory. Upon the introduction of charge, this theory pre- 
dicts the existence of several phases whose structure is 
analogous to the ncmatic and smectic phases of liquid 
crystals^SiiSii^ In the stripe picture, charge and spins 
are spatially segregated with the spins forming antiphase 
domains. In the smectic phase, long-range correlations 
exist producing incommensurate Bragg peaks. In the ne- 
matic phase, the stripes only have orientational order 
and therefore the antiferromagnetic correlations will be 
short-ranged. Also, since the ground state of the stripe 
phase breaks a continuous symmetry there should be no 
spin-gap. This theory seems to give a good account 
of the phase diagram of the La2-2:Sr:cCu04(LSCO or 
214)^ systems which have been studied in great detail. 
For YBCOe.e Mook et alm^ obtained the first evidence 
for stripes along the b axis. The overall picture in the 
YBa2Cu306-i-a; system is still not clear, however, since 
a complete study as a function of doping has not been 
completed and, as we will also show, the structure of the 



2 



8000 



S 6000 

3 
O 

U 

■t 4000 

S 



2000 



(0 6) Bragg Peak 
Rocking Curve 



YBCO Ortho-II 
[12 29 51 120] 




-2.4 



-0.8 0.! 
CO (degrees) 



FIG. 1: Rocking curve of the six-crystal composite sample at 
the (0 6) Bragg reflection. The collimation from source to 
detector is given in minutes of arc within square brackets. 



spin dynamics is highly dependent on structural disorder 
and impurities. 

The most striking feature of the spin spectrum in 
the YBCOg+a system is the presence of an intense 
resonance peak at low temperatures vsrhich is well 
defined in both momentum and energy^ The reso- 
nance peak has been observed in the YBa2Cu306+x, 
BiaSraCaCuaOg+a (BSCCOg+a) (Ref. El), and the sin- 
gle layer TlaBaaCuOe+s (Ref. lUh systems suggesting 
that the resonance is a universal feature of all cuprate 
superconductors. Many theories have been developed to 
explain the resonance peak in terms of the spin-fermion 
and band structure models with the Cu^"*" spins being 
coupled to the quasiparticlesiSfliSiiSS Since these theories 
predict that the resonance is a direct consequence of a 
gap opening in the quasiparticle channel, no well-defined 
resonance is predicted to exist in the normal phase. 

We present a comprehensive study of spin fluctuations 
in the oxygen ordered ortho-II YBCO superconductor in 
both the normal and superconducting states up to ~ 40 
meV energy transfer. The oxygen doping in the system 
studied YBCOe.s corresponds^^ to a hole doping in each 
CuOa plane of p=0.09. We expect that the high degree 
of structural order will remove some of the scattering and 
damping seen in the spin response from samples where 
there is disorder from twinning, poor chain order, and 
short oxygen chain segments randomly located. We be- 
lieve that much of the previous discrepancies in the neu- 
tron scattering data can be reconciled by the existence of 
such disorder and a consistent picture of the undcrdopcd 
region can be formed. 

This paper is divided into parts that deal with the 
normal and superconducting phases. The section on the 
normal phase will present new results on the temperature 
dependence of the low-energy excitations that show that 



no spin gap exists and that one-dimensional incommen- 
surate scattering is present that obeys w/T scaling. In 
the section on the superconducting phase we show how 
the low energy scattering is suppressed, and that a reso- 
nance peak grows that is much better defined in energy 
in YBCO ortho-II than it is in disordered YBCO. Wc end 
with a discussion of the integrated intensities and the to- 
tal moment sum rule, where we determine the absolute 
weight of the resonance in the superconducting phase and 
the surprisingly large fraction of resonance weight that is 
already present in the normal phase. Comparisons are 
made with the pseudogap phenomenon and with theory. 



II. EXPERIMENT 

The sample consisted of six orthorhombioSiiS^ crystals 
of total volume ^^6 cm"^ aligned on a multi-crystal mount 
with a combined rocking curve width of 1.5° as shown 
in Fig. n The lattice constants were measured to be 
a=3.81 A, &=3.86 A, and c=11.67 A. Details of the crys- 
tal growth, the detwinning by stress along the a axis and 
the oxygen order have been given earlier iSSiSi The (2 
0) scan of Fig. |21 shows that the majority domain oc- 
cupies 70% of the sample volume. The peak at larger 
IQI (11=2) is the (2 0) Bragg peak from the majority 
squeezed domain, and the peak 11=1.98 is the (0 2 0) of 
the minority domain. An independent check is obtained 
from the satellites produced by oxygen chain order at 
(3/2 0) and (0 3/2 0). Fig. El shows scans through the 
oxygen superlattice peaks. Fits to resolution-convolved 
Lorcntzians showed that the oxygen correlation length 
exceeded ~100 A in the a and b directions, while it was 
approximately 50 A along the c direction. The correla- 
tion lengths compare favorably to the x-ray characteriza- 
tion of highly ordered ortho-II indicating the high quality 
of our crystals 

The measurements were carried out with the DUAL- 
SPEC triple axis spectrometer at the C5 beam of the 
NRU reactor at Chalk River Laboratories. A focusing 
graphite (002) monochromator and a graphite (002) an- 
alyzer were used. A pyrolytic graphite filter in the scat- 
tered beam eliminated higher order reflections and the 
flnal energy was fixed at Ej = 14.6 meV. For unpolar- 
ized measurements the horizontal collimation was set at 
[33' 29' S 51' 120'] for energy transfer greater than 10 
meV and set to [33' 48' 5" 51' 120'] for energy transfers 
below 10 mcV where the magnetic scattering is weaker. 
The vertical colhmation was kept fixed at [80' 240' S 214' 
430'] . The six-crystal assembly was mounted in a closed- 
cycle refrigerator that was carried on an open C-cradle. 
For constant-energy scans along the in-plane directions 
the rotation axis of the C cradle was along the [001] di- 
rection. This permitted scans along (H H 0) in zones 
(1/2 1/2 L), and, by rotating the C-cradle by ~26°, also 
along (H 3H 0) in the (1/2, 3/2, L) zone. To study 
the incommensurate nature of the magnetic Q-correlated 
peaks, scans were made along the [H 0], [0 K 0], and 
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FIG. 2: Radial scan tiirougli tiie (2 0) Bragg peak. Tlie 
Gaussian fit sliows tliat tlie majority domain occupies 70% of 
tlie sample volume. The magnetization inset shows a sharp 
superconducting transition temperature at 59 K with a width 
of - 2.5 K. 



[0 L] directions independently with respect to the zone 
center (1/2, 3/2, L). To do this, the rotation axis of the 
C-cradle was placed along the [TlO] direction. By slaving 
the C-rotation to the L value we could access a general 
(H K L). 

To confirm the magnetic origin of the scattering we 
made polarized measurements with a Heusler (111) 
monochromator and analyzer. A flipper coil was placed 
on the incident side and a graphite filter was used on the 
scattered side to filter higher order neutrons. A Mezei 
flipper in the incident beam allowed spin-flip (SF) and 
non-spin-flip (NSF) cross sections to be measured. At the 
sample two pairs of coils applied either a weak horizontal 
field (HF) parallel to Q or vertical field (VF) perpen- 
dicular to Q (of — 3-5 gauss) to control the direction of 
the neutron spin at the sample. With the flipper on, the 
difference between the HF and VF count rates gives the 
spin-flip scattering from the magnetic electrons alone in- 
dependent of nuclear incoherent or phonon scattering^ 
To prevent any possible depolarization from the Meissner 
state the sample was always warmed to 100 K (—40 K 
above the transition temperature) and then field cooled 
to low temperatures. During scans the field direction 
with respect to the sample never changed more than a 
few degrees thereby minimizing the effects of depolariza- 
tion on the neutron beam. 

We measured the polarized scattering with two differ- 
ent Heusler analyzers. The first small analyzer, 7.9x3.2 
cm^ , was used for the constant energy scans. The hori- 
zontal coUimation was set to [33' 48' S 126' 120'] and the 
vertical colhmation was [80' 72' S 120' 430']. The second 



FIG. 3: Radial scans through the oxygen ordering superlat- 
tice peaks at (3/2 0) and (0 3/2 0) (units of a*). The 
horizontal bar indicates the resolution. A fit to a resolution- 
convolved Lorentzian shows the oxygen ordering correlation 
length exceeds 100 A. The larger peak (filled squares) is from 
the majority domain and the lesser from the minority domain. 



taller analyzer of dimensions 13x12 cm^ was installed for 
the constant-Q studies through the resonance. The hor- 
izontal coUimation was then set to [33' 48' S 126' 120'] 
and the vertical colhmation was [80' 72' S 280' 430']. The 
hipping ratios from both the small and large analyzer 
crystals were measured at the (2 2 0) Bragg position for 
Fii = Ef = 47.6 meV, the incident energy at the reso- 
nance peak. The flipping ratio with the small analyzer 
was found to be 19:1 for both VF and HF fleld directions. 
The flipping ratio for the tall analyzer was measured to 
be 15:1 for the HF and 12:1 for the VF directions. By 
using a taller analyzer we were able to obtain a factor of 
two increase in the resonance scattering per unit time. 

To compare our results with those of other groups and 
to connect with theory we have put our measurements of 
x"(Q,u;) on an absolute scale. We have calibrated the 
spectrometer by measuring the integrated intensity of an 
acoustic phonon near the strong (0 6) Bragg peak. By 
comparing the measured energy-integrated intensity to 
that calculated in the long- wavelength limit, where eigen- 
vectors are known, we were able to obtain the calibration 
constant that puts the magnetic intensities on an abso- 
lute scale. Details are provided in the appendices. The 
appendices also show how we made the correction for 
monitor contamination by higher-order wavelength neu- 
trons. This correction is particularly important at small 
energies, where its neglect may lead to an underestimate 
of the low energy response. 
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FIG. 4: The peak dynamic susceptibility, x", at Q = (tTiTt) 
as a function of energy transfer for energies below 16 meV 
in the normal and superconducting states. In the normal 
state scattering is present at the lowest energies studied and 
increases roughly linearly with energy transfer (line). In the 
superconducting state the scattering is suppressed but not 
eliminated for energies below 16 meV (the lines are a guides 
to the eye) . 



III. THE NORMAL STATE 



A. Gapless normal state spin excitations 
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FIG. 5; Normal state neutron scattering scans through (1/2 
3/2 2) at energy transfers of 12.4, 8.3, and 6.2 meV. The 
correction for higher order contamination of the monitor rate 
has been made. In the left two panels the scattering along 
[100] reveals incommensurate features, but along [010] it does 
not (see text). The solid curve is a fit of the spin- wave model 
of Chou et al discussed in the text. Their relative weight is 
given by the 70:30 domain ratio. The broken lines give the 
linear background. 



The existence or not of a normal state spin gap has 
many important theoretical implications. Theories in- 
volving a broken discrete symmetry such as orbital cur- 
rents require the existence of a normal state spin-gap 
at low temperature. In one model, the d-density wave 
(DDW) theory, the orbital currents are equivalent to 
static weak moments with Ising-like behavior, and the 
theory predicts that a magnetic Bragg peak will appear 
in the underdoped phase at temperatures well above the 
onset temperature of superconductivity>i& In this model, 
in which the orbital currents are in opposite directions in 
adjacent unit cells, the DDW peak is predicted to break 
the spatial symmetry of the copper oxide plane and to 
appear at Q=(7r,7r). In the fluctuating orbital current 
theory of Wen and LeeSiS there is no static order, and 
there is a crossover to fluctuations that compete with 
superconducting fluctuations in the normal phase. It is 
possible that this might still lead to a suppression of the 
magnetic spectral weight well above Tc- 

Experiments on underdoped YBCO in which no ob- 
servable scattering below ~ 10 meV has been detected 
have been taken as evidence for a pseudogap phase in 
the normal stated In terms of this interpretation the 
pseudogap is defined to be a range in energy transfer 
where there is zero or negligible spectral weight in the 
energy response as measured by the dynamic susceptibil- 
ity, and which lies below a range with appreciable 
weight. Another definition of the pseudogap phase arises 



from NMR studies which have shown the suppression, 
not absence, of both the Knight shift and the relaxation 
rate, l/TiT, well above the superconducting transition 
tcmperaturei2i22iSi24 Based on this we interpret the nor- 
mal phase pseudogap to be a spectral range in which 
cither there is negligible scattering or there is substantial 
suppression of the scattering with decreasing tempera- 
ture. Therefore, in terms of neutron scattering, evidence 
for a pseudogap must come from both the energy and 
temperature dependence of x". Recent work on near op- 
timally doped LSCO has shown some evidence for a spin 
pseudogap based on both the energy and temperature 
dependence of the scattering!^ As regards the observed 
dynamics we find no evidence for a gap or pseudogap in 
the normal phase of YBCOg.s based on the temperature 
and energy dependence of x" ■ In this section we present 
our results on the energy dependence of x" pointing to 
gapless excitations in the normal phase. The temper- 
ature dependence showing x" increasing in the normal 
phase with decreasing temperature is presented in later 
sections. 

In an extensive survey of the low-energy spectra we 
have observed that scattering centered on Q = (tt, tt) 
clearly exists in the normal state and extends to the low- 
est energies studied 4 meV) . A summary of our re- 
sults, shown in Fig. ^ indicates the absence of any clear 
normal state spin pseudogap in contrast to the suggested 
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behavior by previous studies. The data were obtained by 
conducting fits to the Q-correlated peak above the back- 
ground that hes under the magnetic correlations centered 
on (tt, tt). We have removed the temperature-dependent 
Bose factor and the anisotropic Cu^+ form factor from 
the Q-correlated scattering to obtain the dynamic sus- 
ccptibihty as discussed in detail in the appendices. 

While the spectral weight at low energy is lower than 
that at the resonance energy (see later), this, in itself, 
is not sufficient evidence for a gap. We have recently 
discussed this result in the context of orbital currents 
and have argued that the presence of low energy exci- 
tations is evidence against the presence of a static d- 
density wave or orbital current phased From Fig. 01 we 
see that the susceptibility decreases monotonically with 
energy below ^15 meV. The clear suppression of scatter- 
ing in the superconducting state will be discussed later. 
Clearly the difficulty of detecting the normal-phase scat- 
tering increases at low energies as the weight declines, 
but the trend is uniform and gives no sign that the spec- 
tral response drops off in a way that might suggest a 
gap. This linear decrease of susceptibility with energy 
transfer is consistent with over-damped spin waves or 
fermi liquid behavior. It is also consistent with several 
early studies of disordered YBCO in the underdoped re- 
gjQj^i36j37j38_ jsfote that the spin gap referred to in several 
studies2i4^2ii& (see Figs. 9, 26 and 13 respectively) is the 
gap in the superconducting state, not the normal phase 
gap whose absence we address here. 

Some typical scans at low energies along the [100] di- 
rection through (1/2 3/2 2), are shown in Fig. |S1 This 
zone boundary is equivalent to the 2D wave vector (tt, tt) 
with L=2 near a maximum of the bilayer structure fac- 
tor. The results are corrected for the effect of higher or- 
der contamination of the monitor rate. If this correction 
were ignored it would lead to an underestimate of the low- 
energy intensity by a factor ~ 1.5 at 6 meV as discussed 
in the appendices. The solid curves represent the over- 
damped spin-wave model of Chou et alm^ (the following 
section) which gives an excellent fit to the data and also 
predicts that the scattering should increase roughly pro- 
portional to the energy transfer huj. Because it is more 
difficult to detect this signal as the energy decreases, it 
is possible that experiments where the sensitivity and 
statistics are optimized for the intense resonance peak 
might have failed to observe the low-energy response and 
so might have led to the inference that a gap was present. 
All of our results indicate that there is no well defined 
normal state spin-gap in Ortho-II ordered YBa2Cu3 06.5. 



B. Incommensurate dynamic structure in 
YBa2Cu306.5 Ortho-II 

To establish the incommensurate structure we made 
purely H and K scans (see Fig. [3 left hand panels) in 
a configuration with good in-plane resolution, where the 
coarse vertical resolution was only 20-30° from the (0 



L) direction and so did not integrate over the incommen- 
surate peaks. We see that the profile shows incommen- 
surate peaks along a* but only a commensurate feature 
along b*. From the H-K anisotropy of the low-energy 
scans in Fig. [Slwe can show that the dynamic spin fluc- 
tuations are one-dimensional along a* with an incom- 
mensurate wave vector 6 ^ 0.06 r.l.u. To do this we 
recognize that the scattering observed is the sum of that 
from the majority and minority domains weighted in the 
oxygen-order ratio of 70:30. Because of the low weight 
of the incommensurate peaks of the b* minority domain, 
the [010] scan shows only a commensurate ridge as the 
scan passes between the strong a* peaks at ±S where 
the resolution picks up its greatest contribution from the 
wings of the broadened a* peaks. 

Since it gives a good description of the scan profile, 
we fit the over-damped spin-wave model of Chou et alm^ 
and thereby extract (5~0.06 after convoluting with the 
spectrometer resolution function. In the Chou model 
linear branches of damped spin waves rise isotropically 
in the 2D plane from incommensurate wave vectors dis- 
placed by ((5, 0) in r.l.u. from the antiferromagnetic point 
(1/2, 1/2), equivalent to Q=(7r, tt)^. Similar behavior is 
found for low-energy excitations in the stripe phase of the 
nickelatea^2iiM4il^ where the stripes are static instead of 
fiuctuating. As we apply the model to the a* domain, q 
is measured from Q and the incommensurate wave vector 
is (line = ((57 0,0), so that the Chou model is the sum of 
scattering from two wave vector origins: 



^(q,.)|,,„^ - ^["(^) + l](^^T(^±^-F) 
where 

hujq± = 7ic(q±qj^c) (2) 

and 



F = hcK = — (3) 

so 

where c is the spin-wave velocity and can be inter- 
preted as a dynamic correlation length. We emphasize 
that the parameter is not strictly the instantaneous 
correlation length as we are not integrating over energy. 
To account for the minority domain we add to Eq. ^ a 
similar term in which the scattering emanates from wave 
vectors qi„c — (0,6,0). Therefore the total cross-section 
is given as follows: 

S'(q,w) = 5'(q,cj)|q._^^=(5_o,o) + RS{q,u})\c^,^^=(o^6,o)■ (4) 
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The ratio of the two ampHtudcs, i?, from the two domains 
is set to the known domain ratio of 30:70 determined from 
the oxygen chain superlattice peaks. After folding with 
the resolution function the fit to the data is performed. 

The excellent agreement displayed in Fig. O shows 
that the slow spin correlations are incommensurate in Q. 
There is no observable static incommensurate peak. At 
these low energies, phonon contamination is unlikely and 
as a check we have verified that the origin of this scat- 
tering at each energy transfer is truly magnetic, from 
its temperature dependence around Tc (sec future sec- 
tions regarding temperature dependence). By iterating 
through each data set we have chosen the spin-wave ve- 
locity, he ^ 300 meV-A, a dynamic correlation length, 
20 A, and incommensurability, (5=0.06, that applies 
at all energies and temperatures. 

As shown in Fig. |31the fit from this model describes 
the line shape along both the [100] and [010] directions 
very well. The minority domain is not strong enough to 
be seen as a separate peak in the [010] scan (lower left 
panel) , both because of its small weight and because the 
majority peaks are broadened by ^ and resolution so that 
they contribute most of the scattering at the commensu- 
rate point q=0. Thus the [010] scan shows a single peak 
at K=1.5 as the scan crosses a saddle point between the 
two majority peaks at (± (5, 0, 0). 

Our analysis provides strong evidence that the low- 
energy spin dynamics reflect a one-dimensional incom- 
mensurate structure along a*. Its dynamics can be asso- 
ciated with the Doppler shift for a right and left moving 
incommensurate spin density wave. Similar incommen- 
surate scattering has been previously observed for oxy- 
gen concentrations of a; = O.Si& and x = 0.7— at en- 
ergies of 20 meV, where the incommensurate wave 
vector and doping is larger. The one-dimensional in- 
commensurate scattering is thus a common feature of 
the YBCOg+K system in the underdoped region. One- 
dimensional incommensurate scattering has also been ob- 
served in Lai. gsSro.osCuO^i^ indicating that it is possibly 
a general feature of all cuprates. 

Physically, our data can be interpreted in terms of one- 
dimensional antiphase domains parallel to the b* chains. 
This may be consistent with the stripe picture^ in which 
the spin correlations are antiferromagnetic in a domain 
between charged walls separated by £ spins along [100], 
where the carriers concentrate. Across a wall the cou- 
pling is ferromagnetic so that the next spins belong to 
an antiphase spin domain tt out of phase with the first. 
Antiferromagnetic coupling across a row of spin vacancies 
would also give antiphase domains— A nice feature of the 
stripe model is that the local spin correlations within a 
domain remain commensurate, as suggested by analysis 
of NMR datai^ Nevertheless, the maximum in a scatter- 
ing experiment moves away from tt to an incommensu- 
rate wave vector, 2it/2£ i.e., 5=1 /2£, the inverse of the 2£ 
spin repeat distance for long-ranged stripe correlations. 
Our data are consistent with the stripe model where the 
charge domain is of length £=6 spins, and the domain 
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FIG. 6: The agreement of the stripe model (lines), in which 
antiphase domains of length £=6 with 6 = 1/2^ = 0.08 occur 
with a correlation length of 27 A, with the observed scattering 
near the antiferromagnetic wave vector (1/2 3/2 2). 



correlations fall off over a length of 7 spins (~27 A) as 
shown in Fig. This fit was conducted by convolving 
the stripe model with the Q resolution only and taking 
the energy dependence to have a simple Lorentzian form. 
One of the predictions of a sharp domain boundary is the 
presence of higher-order satellite peaks which are clearly 
washed out by the resolution. This model gives a (5 of 
0.08, slightly higher than that deduced using the Chou 
model but consistent given the errors of the fit and the 
fact that we have neglected energy dispersion in fitting 
the antiphase domain model. 

Because our data are consistent with the idea that an 
incommensurate structure exists in the normal state, it 
is possible that the stripe phase is a precursor to super- 
conductivity in the underdoped YBCOg.s system. The 
one-dimensional nature of the scattering is difficult to in- 
terpret as originating from band structure effects— We 
find that for models such as Refs. the anisotropic 
hopping caused by the small 1% a-b splitting should 
produce little difference in the bands and peaks in x" 
along the a* and b* directions. This is confirmed by 
density-functional calculations for the orthorhombic unit 
cell which predict nearly isotropic hopping integrals in 
the bandstructure— It seems more likely that the intrin- 
sic structure of the spins (and carriers) in a plane favors 
stripe order and that the small a-b splitting as well as 
the filled oxygen chains serves to select a* to orient the 
preexisting stripes— Our results are similar to those 
found in the LSCO system. In particular, the incommen- 
surate wave vector in this experiment is consistent with 
those found in the monolayer cuprates with the same hole 
doping— 

A comparison of the incommensurability with doping 
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FIG. 7: The reduction of the superconducting transition tem- 
perature upon reducing the oxygen doping below optimum, 
as it relates to the position of the incommensurate peak in 
cuprate superconductors. The incommensurate wave vector 
along the Cu-O-Cu direction in reciprocal lattice units (r.l.u.) 
such that the wave vector is 2n5/a. 



for bilayer and single layer systems is shown in Fig. As 
can be seen, our data are consistent with values of qinc 
observed in the LSCO system for similar hole doping. 
Even though data in the YBCO system are currently very 
limited, the good agreement to-date in the underdoped 
region points to a possible universal trend between the 
hole doping and the spin stripe width. Following Yamada 
et al&, Balatsky and Bourges^ showed that the overall 
half-width at half maximum of the Q-correlated peak 
below the resonance, whether it shows incommensurate 
modulation or not, is proportional to Tc- Our results 
are in agreement with this relation. This suggests that 
the widths in YBCO in part arise from incommensurate 
modulations and so can be interpreted within a stripe 
picture similar to that in LSCO. 

The fact that S is small, combined with short corre- 
lation lengths and twinning, could explain why many 
studies in this doping region have observed only a single 
commensurate peak at (tt, tt). This point was the main 
assumption of the analysis of Chou et al. who fit a se- 
ries of commensurate inelastic magnetic peaks to a model 
with spin waves originating, with equal weight, from the 
two slightly incommensurate positions of q=(0,5, 0) and 
q={S, 0, 0)i2& It is also clear from our data that in a sam- 
ple with equally populated structural domains, a single 
commensurate peak would be observed because of the 
small incommensurability. 

The incommensurate correlations along a* are a robust 
feature of the behavior and not a result of superconduct- 
ing order in underdoped YBCO. They are present in the 
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FIG. 8: Illustration of in-plane wave vectors accessible in the 
scattering planes (H 3H L) and (H H L). The actual scans 
were chosen to lie in the (H 3H L) plane so that the vertical 
resolution (ellipse) would integrate over the stronger incom- 
mensurate peaks along a*. 



normal phase up to temperatures of at least ~ 120 K, 
and, as we shall see later, they persist in the supercon- 
ducting phase up to ~ 25 meV (Fig. [T^ . 

In scanning tunnelling microscopy Davis, de Lozanne 
and collaborators^ have observed a charge modulation 
along the b* (chain) direction in optimally doped YBCO 
with a 13 A periodicity cell. This would correspond to a 
reciprocal lattice wave vector of 0.3 r.l.u. In underdoped 
YBC06.5 there is no evidence for this charge modulation 
along b* nor of any double-length spin modulation 0.15 
r.l.u., both of which would have been easy to resolve. The 
spin modulation wave vector lies along a*. However, the 
scanning tunnelling microscope images the chains not the 
electronic structure of the planes. 



C. Temperature dependence of x"(Q? and uj/T 
Scaling in the Normal State 

To study the Q-correlated peak as a function of energy 
and temperature, the sample was aligned in the (H 3H 
L) and the (H H L) scattering planes. This configura- 
tion was chosen so that the vertical resolution integrated 
over the intensity from the two slightly out-of-plane in- 
commensurate peaks (see Fig. (S))- Since the width of the 
Q-correlated peaks was found, within experimental error, 
to be independent of temperature, the peak susceptibil- 
ity was then a good measure of the integrated intensity. 
To extract the susceptibility in absolute units we have re- 
moved the anisotropic Cu^+ form factor multiplied by the 
bilayer structure factor as discussed in the appendices. 
The bilayer structure factor for odd-symmetry fluctua- 
tions will be shown later to give a good account of the 
L-dependence of the scattering. 

The temperature dependence of the susceptibility is 
shown in Fig. |^for energies of 12.4 meV and 24.8 meV. 
One of the most striking features is the presence of a clear 
suppression of the spin fluctuations at 12.4 meV below 
the superconducting transition temperature, as will be 
discussed in the next section devoted to the supercon- 
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FIG. 9: The peak susceptibility at 12.4 meV and 24.8 meV 
as a function of temperature. For the normal phase, the solid 
curves are fits to the uj/T scaling analysis of Birgeneau et al. 
At 12.4 meV the normal state follows the scaling analysis. 
A clear suppression of scattering is observed in the supercon- 
ducting state. The scattering at 24.8 meV in both the normal 
and superconducting states continues to grow on cooling and 
no longer follows the nearly constant temperature dependence 
predicted by cj/T scaling. The dashed line in the lower panel 
is a guide to the eye. The local and low-frequency susceptibil- 
ity sensed by the NMR relaxation rate of ^''Cu is suppressed 
on cooling below a temperature T*, while the susceptibility 
at the antiferromagnetic wave vector increases uniformly un- 
til the onset of coherent superconductivity. NMR data is for 
YBC06.64 taken from Timusk and Statt. 
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FIG. 10: A summary of the x" is plotted for all normal state 
data below 16 meV energy transfer as a function of uj/T. The 
temperature dependence of the low energy scattering is very 
well described by uj/T scaling. 



maximum NMR relaxation to occur at much too low 
a temperature ~30 K and cannot explain the observed 
T* ~ 150 K (estimate via Eq. 21 of Tranquada et alm^). 

In early studies of the temperature dependence, x" 
in disordered YBCOg.s was found to follow a sim- 
ple uj/T scaling rclationi^ The same cu/T scaling was 
very successful in describing the temperature depen- 
dence of the integrated susceptibility in underdoped 
La2_a;Sr3;Cu04i^£*S The analysis assumes that the dom- 
inant energy scale is set by the temperature and predicts 
the following temperature and energy dependence of the 
scattering, 



^^arctanf^ 



a2T3 



(5) 



Here the term x"{^tT = 0) represents the limit of x" 
as the temperature goes to zero. A detailed study on 
disordered YBCOe.s showed^^ that the temperature de- 
pendence of all scattering up to 33 meV is described 
well by just the first term with ai set to 0.9: 



ducting state. In the normal state, x" increases with 
decreasing temperature. As previously discussed this is 
further evidence for the absence of a normal state pseudo- 
gap in this energy range. Most surprising is that x" {Q^ ^) 
continues to grow on cooling while the NMR ^'^Cu relax- 
ation 1/TiT (proportional to x"{^)/'^) declines below 
about T* - 150 K,^ The T* from NMR has been as- 
sociated with a pseudogap, but it does not grow signif- 
icantly with underdoping (Tc < 60 K);^ as would be 
expected for a pseudogap that many believe grows with 
underdoping. We cannot exclude a smaller gap than 4 
mcV (Fig. 0}. but such a small gap would predict the 



The inclusion of higher terms in the expansion was found 
not to noticeably improve the fits. This scaling was also 
originally proposed in the context of the marginal Fermi 
liquid theory to explain many of the normal state physical 
properties of the cuprates»S The marginal Fermi liquid 
theory does not give the strong g-space structure seen 
here. 

For ordered YBCOe.s in its ortho-II state we find that 
in the normal phase this scaling form is in excellent ac- 
cord with the spin response at low-energies. Surprisingly 
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the scaling function fits the data with the same coeffi- 
cient ai=0.9 as for the disordered system. The relation 
is plotted in the normal state for 12.4 meV and 24.8 meV 
transfers in Fig. |^and shows excellent agreement at 12.4 
meV but not at 24.8 meV. The results for energy trans- 
fers below 20 meV are summarized in Fig. 1101 The fig- 
ure was obtained by fitting (with ai=0.9) the uj/T scal- 
ing relation at each energy to derive a single parameter 
X"(w,r = 0), from which x"(w, T = 0) follows 

for all temperatures and energies. We conclude that uo/T 
scaling accurately describes the normal state low-energy 
huj and temperature dependence. 

Despite the fact that the scaling analysis works very 
well for low energies, the scaling previously observed and 
predicted suggests that for higher energy transfers the in- 
tegrated susceptibility should be only weakly dependent 
on temperature. This is clearly not the case at 24.8 meV 
(Fig. inj where the response grows strongly on cooling. It 
also does not hold at higher energies where we see x" hi- 
creasing continuously with decreasing temperature. This 
difference in scaling in the normal state shows a clear 
departure of the scattering in YBCO ortho-II from that 
in the disordered systems previously studied which ob- 
served scaling up to at least 33 meV. We note that a 
breakdown of w/T scaling has been observed at very low 
energies in Lai.g8Sro.o2Cu04 and is thought to originate 
from a small gap in the excitations due to the out-of- 
plane anisotropy^ 

A possible explanation for the presence of one- 
dimensional incommensurate scattering, oj/T scaling, 
and the breakdown of the scaling relation at higher en- 
ergies can be found in the stripe model. The stripes con- 
sist of one-dimensional domains tt out of phase with each 
other. As discussed by Zaanen et alm^ for a stripe liquid, 
where there is no long-range order and stripes are fluctu- 
ating, there are two extreme energy scales. The first is at 
high energy where the domain walls appear static, and 
the second is the hydrodynamic limit where the stripes 
are in a fluid state. For the second hydrodynamic (low 
energy) limit the energy scale of the system is set by the 
temperature. These assumptions underlying tu/T scaling 
suggest that the fluctuating stripe model provides a nat- 
ural explanation for the uj/T scaling we clearly observe at 
low energies in the normal state. It is also argued by Zaa- 
nen et al. based on hydrodynamics that x" ^ oj/T which 
can also be derived from the uj /T scaling relation by ex- 
panding the expression to first order. Therefore stripes 
do not predict the presence of a spin pseudogap in the 
normal state but predict scattering at all energies and 
that x" ~ w at low energies. This is exactly the model 
applied here to our low scattering and originally used by 
Chou et al. This is further illustrated in Fig. 0] which 
shows the scattering at very low energies decreasing as 

LO. 

The stripe model also predicts the breakdown oi ui/T 
scaling. At some point the energy will be sufRciently 
large compared with the fluctuation frequencies that the 
fluid domains will appear almost static. In our line-shape 
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FIG. 11: Constant energy scans of the scattering at 8.3 meV 
energy transfer. The incommensurate behavior persist un- 
changed apart from amplitude upon entering in the super- 
conducting state. The background beyond the peaks and the 
fast neutron background scales as [n(ti^) + l]. The solid lines 
are from double Gaussian fits to the data. 



analysis this energy would be given by lDc ~ ?ic/^o = 300 
mcV •A/20 A= 15 meV. Therefore we would expect that 
for energies above ~ 15 meV the internal domain dynam- 
ics will start to play a role in setting the energy scale and 
lead to a breakdown oi uj/T scaling. This is qualitatively 
what is seen in our data - perfect uj/T scaling up to at 
least 16 meV followed by a clear departure from scaling 
at around 25 meV. 

These qualitative ideas can be extended to explain the 
role of disorder. Disorder would decrease the correla- 
tion length ^ and therefore increase the characteristic fre- 
quency LOc ■ This would provide a natural explanation for 
why early studies on highly disordered crystals followed 
scaling up to high energies ~ 33 mcV. These arguments 
provide some support to the idea that the presence of 
fluctuating stripes in a liquid phase can explain many 
qualitative features of the data in the normal state in- 
cluding the temperature dependence and line shape. 

Lo/T scaling can also be interpreted as the result of the 
close proximity of a quantum critical point^ where en- 
ergy and temperature are interchangeable. Such critical 
points have been suggested to exist in the over-doped and 
near optimally doped region of the generic cuprate phase 
diagramiSiSi The fact that we observe scaling in the un- 
derdoped region and that scaling is observed over a broad 
region in the LSCO system suggests that the scaling is 
not the result of a quantum critical point. This assertion 
does depend on the extent and size of the cross-over re- 
gion. To address the issue of quantum criticality clearly 
a more detailed study as a function of hole doping would 
be required in the YBCOe+x system. 
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FIG. 12: Constant energy scans of the scattering at 12.4 meV 
energy transfer. The hne shape does not change upon enter- 
ing in the superconducting state. The background scales as 
[n(a;)-|-l]. The solid lines are from double Gaussian fits to the 
data. 



IV. SUPERCONDUCTING STATE 

A. Suppression of Low Energy (< 16 meV) 
Scattering 

On entering the superconducting state, the scattering 
at energies less than ~ 16 meV is suppressed but not elim- 
inated as shown in Figs. IHI^Illlandll2l Thus there is not 
a complete gap. Moreover, the incommensurate structure 
remains in the superconducting state. As shown in Fig. 
Othc suppression clearly starts at T^. and cannot be taken 
as evidence for a normal phase pscudogap as previously 
discussed. This suppression at the onset of superconduc- 
tivity clearly departs from the uj/T scaling relation, in 
contrast to what was derived from early studies on dis- 
ordered systems for similar oxygen concentrations. Our 
measurements also agree qualitatively with the result of 
early studies on higher oxygen dopings which showed a 
departure from uj /T scaling and a suppression of scatter- 
ing in the superconducting statcj^ This response is con- 
sistent with the behavior in underdoped LSCO, where the 
low energy scattering has been found to be suppressed in 
the superconducting stateS^ with almost complete gap- 
ping at the lowest energyiSi For YBCOg.s we find that 
the energy below which the response is suppressed in the 
superconducting state is consistent with the gap derived 
by other groups in the YBCOg+k system for higher dop- 
ing^ and with the universal curve of Dai et alw^ (see 
their Fig. 26). 

Thus, although there is no evidence for a normal phase 
gap, there is clear evidence of formation of a gap as su- 
perconducting order develops, with the degree of gapping 
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FIG. 13: Constant energy scans of the scattering at 24.8 meV 
energy transfer showing a clear flat-top line shape indicative 
of incommensurate scattering. This scan was performed in 
the (HHL) scattering plane. The solid line is from a double 
Gaussian fit to the data. 



increasing as the energy tends to zero as seen in Fig. 0] 
We do not observe a full gap, and indeed would not ex- 
pect to do so in a superconductor where the pairing gap 
has nodes, regardless of the gap symmetry. 

In contrast to the suppression of the superconducting 
response below ^ 16 meV, for energies greater than 24 
meV the response is enhanced (Fig. |^. It is interesting 
to compare this observation to simple BCS theory which 
predicts a gap at 2A = S.SksTc ~ 18 meV. Levi n^^i^^a nd 
coauthors have calculated x"(Q, w) in the random phase 
approximation (RPA) and found a partial suppression of 
the response in the superconducting phase. There is qual- 
itative agreement with our results (Figs. 01 and 13 . The 
temperature dependence of x" for both the s-wave and 
d-wave cases have been computed in the random-phase 
approximation by Bulut and Scalapin0f2£ Their calcula- 
tion was done assuming a single band with only nearest 
neighbor hopping. The hole doping was ^=0.15, some- 
what above the doping of ^=0.09 that is estimated for our 
ortho-II YBCO system from its T^. For a gap with s-wave 
symmetry the imaginary part of the susceptibility shows 
a strong suppression or gap at all energies up to Tujj ^ 
2ksT'c. The d-wave model shows a suppression for ener- 
gies below Tiuj ~ ksTc and an enhancement for energies 
of k^Tc and above. The d-wave results capture only the 
qualitative trend of our data, in the sense that they show 
suppression at low energies and enhancement at high en- 
ergies. Quantitatively, however, the model fails to predict 
the energy scale at which the crossover from suppression 
to enhancement occurs, for it predicts E '^ksTc, which 
is ~ 5 meV, whereas we observe crossover at 16 meV 
or ^^SksTc- Despite this disagreement the general trend 
of our data is consistent with the formation of a d-wave 
instead of an s-wave gap which would show a suppression 
of the scattering at all energy transfers. 
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As regards the spatial spin correlations we observe that 
they continue to exist with the same incommensurate 
form in the superconducting phase as they did in the 
normal phase, and this despite the suppression of the sus- 
ceptibility by the superconducting order. The dynamic 
stripe structure is robust not only with respect to tem- 
perature, but also with respect to energy, for it is present 
up to ~ 25 meV as illustrated in Fig. ^| 

In the context of the stripe picture previously discussed 
for the normal phase, we observe that the spin suppres- 
sion in the superconducting phase occurs over an energy 
range similar to that of the hydrodynamic region of the 
striped normal phase, i.e., the energy range where we 
found that uj/T scaling describes the normal state tem- 
perature dependence. This suggests that the supercon- 
ducting order acts to suppress the amplitude of the pre- 
cursor stripe fluctuations of the normal phase, but does 
not essentially change their spatial or temporal character. 



B. Resonance at 33 meV 

One of the most dominant features of the scattering in 
the superconducting state, as shown in Fig. 1141 is the 
presence of an intense and well-defined peak in energy 
at 33 meV energy transfer. We have verified that this 
resonance peak is magnetic in origin from three features: 
it exhibits the clear L dependence (Fig. I15|l expected 
for the bilayer structure factor of the two copper oxide 
planes, its intensity decreases with increasing tempera- 
ture as does magnetic scattering (Fig. I18|l , and it appears 
in polarized neutrons in the purely magnetic difference 
channel, as will be discussed later. 

In Fig. El we show that the bilayer structure factor, 
scaled by a single amplitude, gives an excellent account 
of the modulation of the resonance intensity along the 
[001] direction. The scan in Fig. El was done by con- 
ducting a series of one-point scans on and off the peak 
and subtracting. The fact that the data are negative 
at the minimum of the bilayer indicates the presence of 
a sloping background. We have also verified that the 
scattering follows the anisotropic Cu^+ form factor by 
comparing constant energy scans done at the (1/2, 1/2, 
2) position to equivalent scans conducted at higher |Q| 
Brillouin zones. This shows that the resonance is consis- 
tent with scattering from Cu^+ spins in the planes which 
fluctuate in antiphase between the bilayers (also known 
as the acoustic or odd-symmetry mode). 

We find, as have others;^ that the resonance is com- 
mensurate. Scans along (K/3, K, 2), (H, H/3, 2), and 
(H, H, 5) directions geometries have been used in previ- 
ous studies to investigate the incommensurate scattering 
as a function of energy and temperature S We have de- 
termined that the resonance is commensurate by scan- 
ning along the (K/3, K, 2), (H, H/3, 2), and (H, H, 
5) directions. The latter two scans are sensitive to any 
incommensurability in Q along a* and show incommen- 
surate peaks for energies below 24 mcV. However at the 
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FIG. 14: x" for ^-U energies studied is summarized in both 
the normal and superconducting states. Two methods of 
measuring x" s^re plotted, the solid points being obtained 
from constant-energy scans and the open points from a di- 
rect subtraction of the high-temperature background suscep- 
tibility measured at 250 K. A clear resonance appears in the 
superconducting state and also in the data at 85 K, indicating 
that a resonant feature does persist into the normal state. 



resonance we observe only a single commensurate peak, 
an example of which we display in Fig. 1161 Its width 
(FWHM) in Q, after correcting for resolution, is ^0.17 
A~^, implying a correlation length of 12 A or about 3 
cells. 

In the ortho-II ordered sample the susceptibility shows 
a linear region followed by a clear upward curvature to a 
well-defined resonance at 33 meV (Fig. I14|l with a half- 
width of 3.6 meV (including resolution). In contrast, 
when the system has the same i=0.5 oxygen content, 
but is disordered, the linear region of the local suscepti- 
bility leads to a broad maximum at ~ 15 mcV without 
the upward curvature typical of a resonant mode>S That 
broad spectrum extends to low energies and differs sub- 
stantially from the present well-defined resonance of the 
well-ordered ortho-II system. In another a;=0.5 sample 
with Tc=52 K the peak of the resonance occurred at ~ 
25 meV again without upward curvature^ The large dif- 
ferences suggests that the effect of structural disorder in 
the chains is to damp out the resonance. Indeed the reso- 
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FIG. 15: Constant energy scan through the resonance along 
the [001] direction. Due to the presence of optic phonons a 
series of two point scans were done with data being taken on 
and off the peak. The solid line is a fit to the bilayer structure 
factor. The fact that the subtraction gives negative values 
at the minimum of the bilayer structure factor indicates the 
presence of a sloping background. 



nance is broadened in disordered systems to give a tail on 
its high-energy wing, resulting in a more symmetric peak, 
whereas in ordered ortho-II we find a sharp resolution- 
Umited decrease above the resonance. When the disorder 
is small, as in the present ortho-II crystal, the resonance 
becomes quite long-lived (of order 10 periods from the 
raw data and longer if resolution is allowed for). Also, 
this naturally explains why in early studies the resonance 
was first discovered in optimally doped systems. These 
have all the chains nearly full and so no oxygen disorder is 
present (i.e., ortho-I chain ordering), while for lower dop- 
ing without stress detwinning a much broader version of 
the resonance is seeniSiS As for the effect of disorder on 
the width in Q. we can understand why, in the compre- 
hensive study of Dai et alm^, the sample with the most 
developed oxygen order. YBCOg.ei as determined from 
clear ortho-II oxygen peaks, exhibited a much sharper 
resonance in Q than did other samples. 

This line of thought would suggest that the resonance 
should be strongly broadened in the LSCO system be- 
cause Sr enters the lattice in a disordered manner, there- 
fore having the analogous effect to disordered oxygen 
chains in YBCO. To-date there has been no clear sign of 
a resonance in LSCO. However, there is evidence for dou- 
bling of the strength of the local susceptibility, relative to 
that in the antiferromagnetic insulator, at '^20 meV, just 
above the superconducting gap^ This could be taken as 
evidence that the superconducting order pushes the spec- 
tral weight to a region above the gapsSSi If this feature was 
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FIG. 16: Scan through the 33 meV resonance at (3/2 1/2 2) 
along (H II/3 2), a direction that is sensitive to any incom- 
mensurate modulation along a* such as was observed well 
below the resonance. A single commensurate resonance peak 
is observed whose width exceeds the resolution (horizontal 
bar). 



taken as the analog of the YBCO resonance, its spectral 
form is substantially different since it is broad with a tail 
to high energies, but may nonetheless be consistent with 
the effects of disorder. The resonance has been found in 
monolayer compounds indicating that a bilayer structure 
is not necessary for its presence^ 

Fig. El shows constant energy scans at 33 meV for 
several temperatures. These scans show that the Q cor- 
relations at the resonance energy still exist well into the 
normal state. Fig. 1181 shows the growth upon cooling of 
the peak susceptibility at 33 and 31 meV. The spin re- 
sponse exhibits a very clear upward break in slope at the 
onset of superconductivity and is clearly enhanced below 
Tc. What is plotted is the peak in absolute units derived 
from the constant-energy scans of Fig. El 

To establish that in the normal phase the 33 mcV fea- 
ture is resonant in time as well as localized in Q it is 
important to observe that its width is less than its en- 
ergy in a constant-Q energy scan. To remove the effects 
of phonons we have subtracted a background scan taken 
at 250 K. Since the magnetic scattering, as determined 
from constant energy scans, at energies greater than ~ 24 
meV, essentially disappears by 250 K this method of sub- 
traction should give the full peak susceptibility x"(Qj ^)- 
We have also subtracted an estimate for the temperature 
independent background by rotating the analyzer crystal 
five degrees and counting. As can be seen from Fig. 1191 
where results for three of the temperatures studied are 
shown, the method of subtracting high-temperature data 
does a reasonable job at removing the phonons. There is 
a very intense optic phonon at 20 mcV which does not 
quite subtract out and gives an apparent rise in scatter- 
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FIG. 17: Constant energy scans through the resonance peak 
as a function of temperature. The fits are to Gaussian peaks 
with the fixed width for all temperatures. A clear peak is 
seen in the superconducting state and it persists in the normal 
state until it almost unobservable at 175 K. 



ing around this energy. Nevertheless the energy scans 
show a clear resonance below and reveal that in the 
normal state the resonance remains as a weaker but still 
well-defined feature. 

The lines in Fig. ^| are the results of a fit to a 
Lorentzian and a temperature dependent background. 
We fixed the width at all temperatures to obtain a sta- 
ble fit in the normal state. The quality of the data above 
Tc prevents us from determining whether or not the reso- 
nance broadens in the normal state. The fit, though, does 
suggest that the peak of the resonance shifts to higher en- 
ergy (as shown in the inset to Fig. I19|l with decreasing 
temperature, reminiscent of a slight decrease in damping 
with decreasing temperature. This shift in resonance fre- 
quency with decreasing temperature can be explained by 
the spin-fermion model of Morr and Finest In the spin- 
fermion model the resonance position is predicted to be 
inversely proportional to the magnetic correlation length. 
We should emphasize, however, that we do not explicitly 
measure any temperature dependence of the correlation 
length. A similar result has been obtained by Fong et 
alw^ for optimally doped YBCO. 

The energy integrated resonance intensity as a func- 
tion of temperature from constant Q scans is seen from 
Fig. 1201 to show a clear enhancement below Tc. Two 
methods were applied to estimate the energy integral. 
In one (closed circles) we took the area under a fitted 
Lorentzian. The background fitted under the Lorentzian 
removes from the estimated part of the response 
that grows on cooling in the range 25-33 meV (in Fig. 
^1 compare T=65 K with 175 K). Furthermore we know 
from low temperature data that the resonance carries a 
long tail to low energies which is not accounted for by 
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FIG. 18: The peak intensity as obtained from constant en- 
ergy scans of the scattering at 33 and 31 meV is plotted as a 
function of temperature. Both sets of data show a clear en- 
hancement of the intensity at Tc and presence of a resonance 
peak at Q=(7r, vr) well into the normal state. 



a symmetric Lorentzian fit. Since we know the range 
25-33 meV to be free of plionons in our temperature- 
subtracted data (see polarized data later), we adopted a 
second method for the spectral weight in which we nu- 
merically summed the observed response from 25 to 43 
mcV (open circles in Fig. I20() . This method will account, 
as we see, for less than the full integral at low tempera- 
tures. 

The fraction of the full superconducting resonance in- 
tensity that remains in the normal phase just above Tc 
(Fig. 120(1 is 25% from the Lorentzian fits and 70% from 
the numerical integration. Either estimate is a substan- 
tial fraction of the low-T resonance. We conclude that 
the resonant feature exists in the normal phase as a tem- 
poral and spatially correlated feature, suggestive, as dis- 
cussed later, that superconducting fluctuations persist in 
the normal phase. A normal phase resonance at 34 meV 
was observedS^ below 150 K in YBCOg.e and used to 
define a pscudogap temperature. As a fingerprint for 
pairing it appears that of order half the supcrconduct- 
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the Bose factor removed. The rise around 20 meV is due to a 
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in our analysis, giving a sloping background. The data shows 
a clear resonance in the superconducting state. The 65 K data 
indicate that a weakened resonance persists into the normal 
state. The lines are a Lorentzian on a sloping background. 
The inset shows that the resonance energy decreases slightly 
with increasing temperature. 
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FIG. 20: The energy integrated intensity of the resonance 
as measured in constant Q scans. The integrated intensity 
in absolute units was obtained by summing the data from 
25 to 43 meV, and also by fitting a Lorentzian to data like 
Fig. 1191 The latter method gives a resonance integral that 
reaches a maximum of ^ 19 ^% at low temperature. There is 
a clear increase in the growth rate of the resonance intensity 
on entry to the superconducting state. Both methods show 
that a substantial fraction of resonant precursor is already 
present in the normal state. 



ing density, albeit incoherent, has already formed in the 
normal phase. 

By combining both methods, Q-scans and constant- 
Q energy scans using temperature subtraction, we have 
been able to arrive at a complete picture of the magnetic 
spectrum up to 40 meV as shown in Fig. 1141 in both 
the superconducting (8 K) and normal (85 K) states. 
The spectrum is obtained at high-resolution because we 
used a low final neutron energy of 14.6 meV for even the 
largest energy transfers. The absolute calibration was 
determined against the known cross-section of an acous- 
tic phonon assuming a paramagnetic (isotropic in spin) 
cross-section (see appendices). The two methods for de- 
termining x"(Q, w) agree very well, further indicating 
the validity of the assumptions used in the subtraction 
analysis previously discussed. 

In Fig. I^we plot the dependence of the low tempera- 
ture resonance frequency on doping (derivecffi from Tc) 
as determined by different groups. It can be seen that 
the resonance frequency in our sample is consistent with 
the overall trend. A linear scaling of the resonance en- 
ergy with Tc (dashed line) is consistent with most of the 
data. However the resonance trend also lies close to the 
predictions of the SO (5)2^ and quantum critical theory of 
Chubuko\iS for which the solid line is a fit to ujres oc (p- 
Pc)^/^ with pc=0.055, as measured carefully in the LSCO 
system^. Because its spectral weight has a strong onset 
at Tc for optimally doped, and an accelerated growth rate 
at Tc for undcrdopcd, there is little doubt that the reso- 



nance peak is intimately related to superconductivity. It 
is not yet clear which effect drives the other; the pairing 
could involve the entire very high energy spin- wave band 
while the resonance is a consequent low-energy concen- 
tration of spin response caused by the pairing. Further 
studies, particularly at lower doping, are required to ad- 
dress this point. We will discuss this topic later in the 
context of the total integrated intensity and the sum rule. 

At low temperatures we clearly see the presence of 
both a commensurate resonance peak, and incommen- 
surate scattering at low energy transfers. This coexis- 
tence points to a possible common origin of both the res- 
onance and incommensurate scattering. This topic has 
been discussed by Batista et al. who have analyzed pos- 
sible common magnetic origins of both the resonance and 
incommensurate scattering iSMi In the analysis of Batista 
et al. the incommensurate features are interpreted as 
spin-waves originating from static incommensurate wave 
vectors. A resonance then results when the spin- wave 
branches meet at the (tt, tt) position. This idea connects 
with our data as our line shape found in the unpolarizcd 
data has a clear tail at low energies and is suggestive of an 
umbrella type dispersion where the excitations meet at 
the resonance position. Such a dispersion in the context 
of the incommensurate to commensurate behavior has 
been discussed in studies of optimally doped YBCOi^i 
However, these ideas are speculative for YBCOg+x as the 
incommensurate scattering is broad and we cannot re- 
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FIG. 21: A summary plot of the resonance energy as a func- 
tion of hole doping ■p for the YBCOe+i system. The solid 
line is a fit Wres oc (p-pc)^''^ as predicted by SO(5) theory. 
A similar doping dependence for the resonance energy was 
also predicted by Chubukov et al. based on the presence of a 
quantum critical point. The dashed line is a linear fit. Both 
fits give reasonable descriptions of the data. 



solve the spin-wave branches suggested by this analysis. 
Experiments on stripe-ordered systems, however, give a 
spectrum with a sharp cut off for Q=7r at a maximum 
energjii^, similar in form to the asymmetric spectrum we 
observe for the YBCO resonance. 

It is not clear how the low-energy fluctuating stripes 
relate to the resonance peak. In the normal state wc see 
formation of incommensurate fluctuations below ^^120 K 
or 2Tc, comparable to the temperature below which the 
resonance first becomes clearly observable. On cooling, 
the growth of resonance intensity accelerates on passing 
into the superconducting phase, concomitant with the 
onset of suppression of the stripe fluctuations. Whether 
one or the other is the progenitor remains an open issue. 



C. Polarization Analysis 

To determine if there is any preferred orientation for 
the spin fluctuations we studied the resonance and low- 
energy scattering in three different zones and scatter- 
ing planes. These included the (H, H, L), (H, 3H, L), 
and (3K, K, L) zones. With unpolarizcd neutrons the 
low-energy incommensurate scattering and the resonance 
peak do not show measurable asymmetry between the (H, 
3H, L) and (3H, H, L) directions. This rules out any en- 
hanced spin fluctuations oriented along either the a* or 
b* directions. A comparison of the intensity in the (H, H, 
L) plane with that in other scattering planes suggests no 
preferred spin fluctuations along the c* direction. This 



is the same isotropic spin polarization behavior as in a 
paramagnetic phase. 

Further evidence for this finding comes from polariza- 
tion analysis. For scattering from a paramagnetic mate- 
rial with a horizontal field (HF) along Q, all of the mag- 
netic scattering will appear in the spin flip channel since 
the neutron detects only the perpendicular component of 
the magnetization. For a vertical field (VF) perpendicu- 
lar to Q the magnetic scattering will contribute equally 
to the spin-flip and non-spin-flip channels. If the mag- 
netic scattering is entirely paramagnetic, the scattering 
in the HF channel will then be exactly twice that in the 
VF channel. Fig. shows polarized scans through the 
resonance at three different temperatures. The HF in- 
tensity is always twice that in the VF channel to within 
error, showing that the scattering arises only from the 
magnetic spins and that the resonance is isotropically 
polarized. The fact that in this paramagnetic system 
we see no preferred polarization in the spin fluctuations 
can be taken to be consistent with transitions from a 
singlet superconducting ground state to a triplet excited 
state. We note that a triplet particle-particle resonance, 
coupled by the superconducting order parameter into the 
particle- hole spin susceptibility, is predicted by the SO (5) 
theory of Demler and Zhang — . We note that a triplet 
resonance was an early prediction of Chubukov et al& 

The temperature dependence of the resonance inten- 
sity agrees with that obtained from the unpolarized con- 
stant energy scans. The agreement proves that magnetic 
resonance scattering persists in the normal state at 65 
K. Since the scattering disappears at high temperatures, 
this unambiguously demonstrates the absence of any con- 
tamination from accidental Bragg or phonon scattering. 
This type of contamination has been a problem in previ- 
ous studies of the resonance modeiii^^ 

Fig. 1211 also shows that the background around the 
resonance peak has the same intensity in the HF and 
VF channels. This indicates that there is no measurable 
continuum or diffuse magnetic scattering other than the 
antiferromagnetic correlations centered on the (tt, tt) po- 
sition. This point is further highlighted by the dotted line 
in Fig. 1221 which gives the fast neutron count with the 
analyzer rotated five degrees. Even though this method 
does not obtain the full background, it does put a lower 
bound on the amount of background scattering. From 
this one can see that most of the background is taken 
into account by the scattering with the analyzer turned. 
The rest of the spin-flip intensity in the wings is close to 
that expected from feed-through from the non spin-flip 
channel. This contribution to the scattering in the spin- 
flip channel is equal to the scattering in the non-spin-flip 
(NSF) channel divided by the flipping ratio. This is equal 
to approximately 10 counts for the monitor used in Fig. 
1221 The NSF feed-through and fast neutron rate account 
for almost the entire background. This shows that there 
is no measurable magnetic diffuse or continuum scatter- 
ing around the correlated peak at (tt, tt). Finally, we note 
that in this energy range there is no observable scatter- 
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FIG. 22: Polarized constant energy scans through the reso- 
nance peak at temperatures below and above Tc- The de- 
crease with temperature of the resonance intensity at 8 K 
agrees with that of the unpolarized data, eliminating any pos- 
sibility of phonon contamination. The open circles are taken 
with guide field parallel to Q (HF) and the triangles with 
field perpendicular to Q(VF). The peak intensity in the VF 
channel is roughly half that in the HF channel, which indi- 
cates no preferred polarization of this excitation. The dotted 
line is an estimate of the fast neutron background. The feed- 
through from the NSF channel can be estimated from the 
hipping ratio to be on average ~ 10 counts. When added to 
the fast neutron background, this almost completely accounts 
for the background around the correlated peak. The temper- 
ature dependent high point at H=0.35 also appears in the 
NSF channel and is likely due to a phonon. The solid lines 
describe Gaussian fits to the data. 
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FIG. 23: Polarized measurement of the spectral form of res- 
onance peak obtained by subtracting the vertical field (VF) 
from the horizontal field (HF) data . The 33 meV peak posi- 
tion, the resolution limited high-energy cut off, and the low- 
energy tail connecting to the incommensurate response, con- 
firm the spectral form got from the unpolarized data (dotted 
line is smoothed unpolarized data scaled by a constant fac- 
tor). 



ing from hydrogen. This would have occurred as nuclear 
spin-flip scattering equally in VF and HF channels. This 
attests to the careful mounting of the sample in dry con- 
ditionsS and to the minimal amount of shielded adhesive 
used to seal the sample cans. 

Since it is established that the scattering is param- 
agnetic, we have have been able to independently verify 
that the low temperature spectral line shape derived from 
unpolarized data (Fig. [T^ is correct. We subtracted the 
spin-flip scattering at (1/2, 3/2, 2.2) in HF and VF chan- 
nels and so obtained the magnetic scattering alone, free 
from the phonon contamination which we had to remove 
in analyzing the unpolarized data. The resulting spectral 
form of the resonance shown in Fig. 1231 is in excellent 
agreement with the unpolarized results (the dotted line 
shows the smoothed unpolarized data). In particular, the 
polarized scans confirm the interesting asymmetric line 
shape consisting of a low-energy tail leading to a sharp 
cut-off beyond the peak of the resonance as discussed 
above. 



D. Integral of S(Q,a;) and the Sum Rule 

From Fig. we have estimated the total integrated 
spectral weight of the resonance at 8 K. To do this a 
symmetric Lorentzian profile centered at 33 meV was fit 
to the low temperature energy scan in Fig. 1141 to obtain 
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a half-width r=3.6 meV and a peak intensity of 1400 
/i^/eV. By assuming a Gaussian hne shape along both 
the [100] and [010] directions, a perfect column of scat- 
tering along the [001] direction, and ignoring the bilayer 
structure factor, we estimate the integrated scattering of 
the resonance to be 



Integral of Resonance: YBCO 



TT"^ I dw I <i3g[n(a;) + l]x"(q,w) = 

0.052^1 



(7) 



per formula unit where d'^q is in reciprocal lattice units. 
Our integrated intensity for the resonance compares very 
well to the results of previous studies for similar oxygen 
concentrations. Because of resolution effects, the true 
comparison between different groups and experiments 
comes from the integral of x" iii both momentum and 
energy. Fong et alm^ calibrated their measurements us- 
ing an acoustic phonon and have obtained an integral of 
0.022 fi% (0.069 ^J.%/^:) for YBCOe.s. As noted in the 
appendices, the definition of Ires used in that work is tt 
times larger than the one used here. Dai et alS^ have 
integrated the resonance to find ~ 0.06 /i^ for YBCOe.e- 
This calibration was conducted using the isotropic Cu^+ 
form factor which, given the Q position in that experi- 
ment, would introduce a small ~ 5 % difference between 
our measurements. A summary of the absolute measure- 
ments of the resonance as a function of Tc is plotted in 
Fig. 131 The ILL and lower Chalk River data points are 
the growth of the resonance integral that occurs below 
Tc- The ORNL and the upper Chalk River data points 
are the total resonance integral. 

An increase in the spectral weight under the resonance 
peak with decreasing doping is predicted by the SO (5) 
theory of Demler and ZhangiSS Their theory relates the 
spectral weight in the superconducting state to the dif- 
ference in exchange energies between the normal and su- 
perconducting states extrapolated to T=0, taken to be 
the condensation energy, which is expected to decrease 
with lower doping. As noted by Demler and Zhang this 
is consistent with the electronic specific heat measure- 
ments of Loram et al. who have found a decrease in the 
condensation energy with decreased dopingiS The con- 
densation energy is proportional to the difference in the 
integral of the dynamic susceptibility. For the spin re- 
sponse below 42 meV we find that in YBa2Cu306.5 the 
integral / duj x"(qi ^) (iio Bose factor) is indeed larger 
in the superconducting phase, since it amounts to 17 ± 2 
/i^ at 8K and only 13 ± 2 at 85K. but the difference 
lies within error. 

An increase with reduced doping of the resonance in- 
tensity also arises in the spin-fermion modeli^ With this 
model, and previous normal state inelastic neutron re- 
sults, Morr and Pines were able to predict quantitatively 
a factor of ~2 increase in the integrated intensity from 
optimal doping to YBCOg.s. The results of Fong et al^ 
show that any growth is by a smaller factor. However, 
Fong's integral includes only the component of the in- 
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lower data are the growth in resonance weight below Tc for 
ILL and below 85 K for Chalk River. 



tegrated intensity which develops below Tc- We find the 
same integral if we remove from our data all but the reso- 
nance growth below 85 K (lower Chalk River point in Fig 
I24|l . However, our results for the total resonance weight 
of the normal plus superconducting state, combined with 
the comparable results of Dai et at, show that there is 
a clear growth of total spectral weight with decreasing 
doping. These results lend support to theories, such as 
SO (5) and the spin-fermion model, which predict such 
an increase. 

Ignoring the effect of the chains, and therefore taking 
only two Cu^^ atoms per formula unit, the total spec- 
tral weight integrated over all all energy and momentum 
should be given by the total moment sum rule: 



J duJ J d^Q Siq,Lj) = 2 X ^5" (S+l)g^, 



(8) 



where j d^q is the momentum space integral over the 
correlated peak. In this equation the factor of 2 comes 
from the fact that we are taking two Cu^"*" ions per for- 
mula unit. Here we have assumed that 5(q, has been 
corrected for the bilayer structure factor and the Cu^+ 
form factor (see appendices) . We do not expect the total 
moment sum rule to be strictly obeyed in the supercon- 
ducting cuprates as the holes will destroy a fraction of 
the Cu^+ moments. Given that the hole doping in our 
sample is p ^ 0.1, this effect should be small and we 
expect the total measured moment to be similar to that 
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predicted by the sum rule. Equating the above integral to 
the cross-section for paramagnetic scattering discussed in 
the appendices, where the relation between ^(q, lo) and 
x"(q, oj) is defined, we obtain the total moment sum rule 
for localized spins as the integral over ^(q, a;): 

I = 7r-^ JdujJ d^'qinito) + l]x"(q, lo) = (9) 

For this gives a total integral of 1= 2fi^. The total 
resonance at low temperatures therefore makes up only 
^ 2.6 % of the spectral weight predicted by the total 
moment sum rule. This agrees very well with estimates 
found in other studiesi&S on YBCOg+x and BSCCOg+x- 
The absolute intensity of the resonance has also been 
computed from band structure theories and is in reason- 
ably good agreement with our experiment Despite the 
fact that the resonance is strongly affected by the ap- 
proach and onset of superconductivity, it contributes a 
tiny fraction of the total moment. Its weakness has led 
to the suggestion that the presence of a resonance is not 
crucial to the onset of superconductivity and that it is 
unlikely the resonance can be associated with the medi- 
ating boson of superconductivityi2i24 This of course does 
not exclude the entire high-energy ^^250 meV spin-wave 
spectrum from providing the boson, while the resonance 
is just a redistribution of a tiny fraction of the low-energy 
spin spectrum in response to the opening of a pairing gap, 
an idea that we now examine. 

To understand the origin of the resonance spectral 
weight we compare in the superconducting and normal 
phases the integral from to 42 meV of the two spectra 
of Fig. El multiplied by the Bose factor. The integral, 
which is equal to tt / ^^^(q, uj), was carried out by poly- 
nomial interpolation of the data. If the resonance growth 
has come from energies larger than our 42 mcV window, 
or from other momenta than than AF region, then we 
would calculate an increase in the superconducting phase. 
For 8 K we obtained a / = / dLux"iq, lo) [n{uj) -|- 1] = 17 
±2 fig and we find 14 ± 2 for 85 K. Although there 
is some growth, within error the total weight of the spin 
correlations is conserved. This indicates that the bulk of 
the spectral weight gathering at the resonance energy at 
low temperatures is transferred from from the low-energy 
incommensurate scattering as this becomes suppressed 
in the superconducting state. Since the Bose factor in- 
creases as ~ w"^ for small energy transfers, a natural 
concern with this type of calculation is the weight car- 
ried by x" at very low energies, a region which is difficult 
to characterize experimentally. To assess its influence we 
note that that integral from to 5 meV at 85 K amounts 
to ~ 0.5 representing only 4% of the total weight. 
Thus the the low-energy limit of x" is not contributing 
a large amount and so validates the accuracy of the inte- 
gral and our finding of a low-energy to resonance transfer 
of spectral weight. 



We now extend our calculations to determine what 
fraction of the total moment resides in the entire spec- 
trum, not just the part under the resonance. By in- 
tegrating numerically up to 42 meV, and with the 
same q-dcpendcncc (Gaussian in the a-b plane and con- 
stant along the c direction) we find a total integral of 
TT-i/ duij d^q [n{uj) + l]x"(q, w) ~ 0.06 fi%. We con- 
clude that the spin response below 42 meV carries about 
3% of the spectral weight available from the total mo- 
ment sum rule. Thus the spin response below 42 meV 
captures only a small fraction of the total moment and 
illustrates the importance of studying higher energy spin 
fiuctuations. 



V. DISCUSSION 

We first discuss the contribution to the absolute total 
moment squared of the entire spectrum up to 43 meV. 
We have shown above that this amounts to ~ 0.06 fi^ 
per formula unit. For a single Cu^^ spin this gives ~ 
0.03 fi%, comparable that observed to develop (on cool- 
ing in zero field) (Ref. Isfil) in underdoped Lai_i,Sri,Cu04 
with x=0.1, a doping similar to that for the present 
YBCOg.s. This is one of the many similarities with the 
LSCO system. We note that both LSCO and YBCO 
ortho-II exhibit scattering that extends to the lowest en- 
ergies in the normal phase with incommensurate features 
characteristic of moving stripes. There is a shorter cor- 
relation length in YBCO ortho-II, yet we have shown 
that the incommensurate wave vector tracks the doping 
equally in both systems (Fig. [71). For LSCO, if we asso- 
ciate the resonance with the enhanced response at ener- 
gies above the region where superconductivity suppresses 
x"(Q, uj), then qualitatively there is much similar behav- 
ior in LSCO and in YBCO. 

On cooling below Tc the strong enhancement of spec- 
tral weight in the resonance energy region is predicted by 
RPA calculations of the susceptibility of a single quasi- 
particle band based on BCS theoryi^S*^ The presence of 
resonance intensity at temperatures above Tc is not. We 
believe, as discussed below, that this reflects local super- 
conducting fluctuations above Tc in agreement with the 
evidence for vortices seen through the Nernst effect 

One aspect of the resonance observed here at high 
resolution appears different from the results of previ- 
ous studies of YBCO. Its spectral shape is asymmetric 
with a sharp drop-off above the peak, and with a much 
slower decrease along a wing below the resonance. This 
wing connects smoothly with the previously noted linear, 
low-energy response part of ^"(Qi lo). Several calcula- 
tions ^ predict the opposite asymmetry as does the t- J 
model22i2S. Norman^ has noted that the calculation for 
x"(Q, uj) drops off too fast below the resonance. On the 
other hand the spin-fermion model of Morr and PinesS^ 
gives a spectral shape for the resonance with a tail to low 
energies as we observe. 

The physical mechanism underlying the resonance 
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peak still remains unknown and is a matter of much de- 
bate. Band structure calculations using a single band as- 
sociate the resonance peak with the opening of a super- 
conducting d-wave gap in the excitation spectrum^2i££ 
Since these theories are mean-field they are not able to 
explain the presence of the resonance in the normal phase 
at a similar energy to that in the superconducting state. 
Other theories such as that of Wen and Lee, where d- 
symmetry normal-phase fluctuations are dynamic^ may 
be more appropriate. The QED3 theory, which takes into 
account vortex-antivortex fluctuations within the pseu- 
dogap region, may also account for the resonance we see 
above — 

A heuristic approach to understanding the spectral 
form of the resonance is to interpret it as the result of a 
sampling of the density of states of the d-superconducting 
(dSC) gap function. In this picture the spectral weight 
of the overdamped low-lying spin response is pushed up 
to energies above a varying gap function. This accords 
with our observation of a transfer from low energies to 
the resonance energy as superconducting order develops. 
We find that a uniform sampling of a dSC form of gap 
fimction, Ak=(Ao/2)(cos(fc2;)-cos(A;j^)) gives a good ac- 
count of the initial rise in x" i ^-nd also the cusp-like 
cut off we observe at the maximum energy of the res- 
onanceiSi This might suggest that the spins couple to an 
incoherent charge spectrum so that the momentum con- 
servation between the spin response and charge response 
of a band model is compromised. It would imply that 
the cutoff beyond the resonance might be a measure of 
twice the maximum gap, 2Ao ~ 33 mcV in a two par- 
ticle picture. However, the recent theory of Eschrig and 
NormanS^, which models the ARPES charge spectrum 
of Bi2Sr2CaCu208+5, predicts for p=Q.l that the charge 
pairing gap 2/S.m = 120 meV ~ SOk^Tc is much larger 
than our observed spin resonance energy of ~ Gfe^Tc, 
and that it increases while the resonance energy declines 
with reduced doping. This much larger energy scale pre- 
sumably provides the less metallic environment in which 
a spin resonance can survive with relatively little damp- 
ing. An attractive feature of the model)2S however, is 
that it predicts an increasing degree of incoherence in the 
charged quasiparticlc properties as doping is reduced, a 
finding that lends support to the momentum averaging 
that underlies our heuristic model. 

We have looked for behavior that might account for the 
decline on cooling of the NMR relaxation rate 1/TiT, 
which has been associated with a pseudogap tempera- 
ture of ~150 K. We note that the low-energy dynamic 
susceptibility in the mcV range does the opposite - it in- 
creases on cooling from room temperature right through 
the pseudogap range, and continues to rise until inter- 
rupted by superconducting pairing below 59 K. Since 
the normal response below 16 meV is linear in uj, it fol- 
lows that its slope x"(Q: ^)/^: which amounts to 0.014 
/i^/mcV^ and which is proportional to the NMR relax- 
ation rate 1/TiT, also increases while the NMR relax- 
ation decreases. The different temperature dependence 



of the NMR and neutron spin response, suggests that 
either the extrapolation from the THz neutron to MHz 
NMR frequency range is far from linear, or that it is at 
the other wave vectors sampled by the local NMR probe 
that the pseudogap suppression takes place. 

The only feature in the neutron spin response which 
has an onset near the presumed pseudogap temperature 
is the resonance peak which grows in below 130 K. The 
increasing weight in the resonance may be interpreted 
in terms of a shift of spectral weight away from very 
low energies. In comparing the results in the supercon- 
ducting state with those in the normal phase, we have 
noted earlier that the rapid increase of the resonance 
intensity below Tc can be accounted for by the loss of 
low-energy intensity below huj ~ 20 meV. Such a picture 
cannot account for the normal state, where the neutron 
spin susceptibility inexorably rises on cooling and gath- 
ers extra weight at the resonance frequency well above 
Tc accompanied by growth, rather than suppression at 
low energies. It is possible that the spectral weight under 
the normal state resonance may be acquired from higher 
energies or momenta, well removed from (tt,tt). How- 
ever, since the resonance energy scale is already much 
larger than thermal energies, it is unlikely that there are 
substantial changes in the higher energy spectrum near 
(tTjTt). On the other hand, the depression of the local 
NMR spin susceptibility (Fig. O, which reflects an aver- 
age over momenta, suggests that much of the resonance 
weight is being transferred from low momentum fluctua- 
tions. 

While the spin correlations at all energies remain vis- 
ible to the highest temperatures as a correlated peak at 
(7r,7r), the formation of the stripe modulation only be- 
comes observable below ~ 120 K. This is in a broadly 
similar temperature range to where we sec the resonance 
begin to grow in. The emergence of resonance intensity 
and stripe modulations are associated in the model of 
Batista et alm^ We emphasize that Batista's theory is for 
an incommensurate static ordered structure, which is far 
from the present situation. 

Calculations for the 2D Hubbard cluster modelS^ have 
led to the conjecture that uniform rather than antifer- 
romagnetic spin fluctuations are causing the pseudogap. 
However, they find no evidence for pairing fluctuations 
for Tc<T<T*, and reject the notion of preformed pairs in 
favor of RVB spin-charge separation. Many of the prob- 
lems of the single-band theories may indeed be overcome 
by a theory (such as spin-charge separation) where the 
dynamics are not described in terms of individual quasi- 
particles^ Spin-charge separation is also favored in the 
theory of P. Leem. 

In the energy spectrum the pseudogap is clearly seen 
as noted above in the charge response given by the opti- 
cal conductivityiSiS^ For YBCOg.T the c-axis conductivity 
is suppressed below an energy of ^ 50 meV. It grows in 
at high temperatures ~ 3Tc. For YBCOg.e the suppres- 
sion of the a-b optical conductivity is below 80 meVi^ 
The connection between the resonance in the spin sus- 
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ccptibility and the charge dynamics has been made in 
the work of Schachinger et al. and Carbotte et al.m^^ 

The resonance has been described as a two-particle 
stated or as a spin excitoni2& Demler and ZhangS^ have 
described in the SO (5) theory how a particle-particle (p- 
p) triplet resonance arises from paired carriers. The spin 
resonance at 33 meV then gives the maximum energy, 
2Ao to create a pair when the d-wave superconducting 
(dSC) gap function is A=(Ao/2)(cos(k2:)-cos(ky)). The 
pair resonance can only be seen in neutron scattering 
when it is coupled to the spin susceptibility by the su- 
perconducting order parameter and hence not above Tc. 
A reduction of the superconducting order by a magnetic 
field would then lessen the coupling to the spin chan- 
nel as found by Dai et al, who observed the resonance 
intensity to weaken with fieldiS^ However Tchernyshyov 
et alm^ have shown that for flat bands near (tt, 0) the 
neutron resonance is more likely a spin exciton near the 
lower energy of Aq. The triplet symmetry of the spin 
resonance was predicted by the quantum critical theory 
of Chubukov et alm^, and also in SO (5), and is confirmed 
by our polarized neutron results which can be interpreted 
in terms of a singlet ground state with a triplet excited 
state. 

It is more difficult to understand the presence of the 
spin resonance in the normal phase, albeit weakened but 
of similar spectral form. Indeed, most theories predict 
that the resonance should be seen in the superconducting 
but not in the normal phase above Tc. Our observation of 
the resonance above Tc does not support these models. 
The quantum critical spin-fermion model of Abanov et 
almSS, (Eq. 140 with a;s/=33 meV), including hot spots, 
does give a broad peak in the normal phase spin response. 
It is a relatively featureless Lorcntzian spectrum, how- 
ever, falling slowly above a spin-fluctuation energy. The 
model does not predict the more peaked spectrum with a 
faster high-energy fall off near the resonance energy that 
we see in the normal phase (Fig. 1141 lower panel). The 
spectral shape of the normal-phase spectrum appears 
more like a weakened version of the superconducting res- 
onance spectrum and provides circumstantial evidence 
that local superconducting fluctuations exist above Tc. 
Thus the neutron response supports the evidence that 
fluctuating vortices exist in a YBCOg.s sample (Tc=50 
K) as has been obtained from the Nernst effect It is 
possible, we suggest, that the normal phase fluctuations 
may be of the correct symmetry to let a local resonance 
appear in the neutron spin response. Since the neutron 
spectrum is independent of the sign of Aq, fluctuations 
of the gap parameter on a timescale slower than the in- 
verse resonance energy should allow a resonance feature 
to remain in the normal state. The density of states 
of the dSC gap function is strongly peaked at its maxi- 
mum energy, and this characteristic energy can persist in 
the normal phase. The loss of pair coherence on heating 
through the dSC superconducting transition can occur in 
the nodal regions, without eliminating the high-energy 
resonance spectrum. By contrast in an s-wave supercon- 



ductor the nodeless order parameter is non-zero every- 
where, and the gap function would have to collapse in 
all directions to effect the transition. Extensions to the 
theoretical models to treat fluctuations so as to extend 
their validity to the normal phase would be valuable. 



VI. CONCLUSION 

We have shown that the antiferromagnetic spin fluc- 
tuations in YBCOg.s ortho-II are gaplcss, with incom- 
mensurate modulation at low energy consistent with dy- 
namic stripes. The spin susceptibility curves upward 
with energy to a well-deflned but asymmetric resonance 
peak at 33 meV followed by a precipitous cut-off. A 
weakened image of the resonance is readily observed in 
the normal phase as a well-defined temporal and spatial 
structure with an asymmetric spectrum similar to that 
of the superconducting state. On cooling, the resonance 
at 33 mcV (~^6.5 k^Tc) and low-energy spin fluctuations 
strengthen until the superconducting phase at 59 K is 
reached. As coherent superconducting order is estab- 
lished the resonance accelerates its growth rate, while 
the low-energy fluctuations below Tc, are sharply sup- 
pressed but not eliminated as expected for a d-wave gap, 
and transfer much of their weight to the resonance re- 
gion. The spin susceptibility, which is confined to nearly 
antiferromagnetic momenta, shows no evidence for an en- 
ergy gap in the normal phase, and we suggest that exper- 
iments in which a pseudogap suppression is seen sample 
a range of lower momenta. However we find that the 
resonance and incommensurate response make their ap- 
pearance below a temperature as high as T*~ 2Tc com- 
parable with the pseudogap temperature. The fraction of 
the low-temperature resonance weight that develops be- 
fore entering the superconducting phase is surprisingly 
large. We suggest that this indicates that substantial su- 
perconducting pairing fluctuations occur in the normal 
phase of underdoped YBCOe.s. 
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VII. APPENDIX A: ABSOLUTE UNITS AND 
CALIBRATION OF SPECTROMETER 

To compare our data with that of other groups and to 
theory we have put our measurements of x"(Q,Li;) on an 
absolute scale by comparing the measured magnetic in- 
tensity to the integrated intensity of a transverse acoustic 
phonon near the (0 6) Bragg peak, e.g. with Q ==(0.15, 
0.15, 6). The use of a phonon from the sample provides a 
good internal calibration which is independent of errors 
resulting from impurities or domains which may be an 
issue with the use of a vanadium standard. 

For fixed final energy and counting time determined 
by the fixed counts in an incident beam monitor whose 
efficiency varied with incident wave vector as 1/ko, the 
measured intensity, I, is directly proportional to the mag- 
netic or phonon cross-section, 



/(Q,w) oc S{Q,uj) 



(10) 



The constant of proportionality (A) can be determined 
from the measured integrated intensity, /(Q)= J duj I(Q, 
uj), of an acoustic phonon. In the long wavelength limit 



h 



[l + nK)]|F^.pQl^e- 



2W 



,(11) 
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where M is the known mass of the unit cell, e" 
is the Debye- Waller factor, [1 -I- n{ijS)\ is the Bose factor, 
i^ATp is the static structure factor of the Bragg reflection 
nearest to where the acoustic phonon is measured, and 
/3 is the angle between Q and the phonon eigenvector. 
Inserting the measured phonon frequency allows a direct 
measurement of the calibration factor A and therefore 
allows any measured intensity to be put on an absolute 
scale. 

For the magnetic scattering we relate the measured 
scattering to the magnetic correlation function bjiiSi 



afi 

5'Q/3(Q,t^). 

The correlation function is related by the fluctuation dis- 
sipation theorem to the imaginary part of the spin sus- 
ceptibility x(Q, w) by 



5Q/3(Q,w) = 7r [ri(a')4-lj ^-^ 



(13) 

Since the paramagnetic scattering is isotropic in spin, 
X" = = Xyy = we can extract x"(Q,w) from 
the measured intensity, I, through the equation 



/(Q,c.)=A^^/2(Q)i?2(Q^)e- 



2W 



^^^(2X"(Q,-)). (14) 
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FIG. 25: The fraction of neutrons hitting the diffracted beam 
monitor is plotted as a function of energy transfer for A/n = 
A/1, A/2, and A/3. The final energy is taken to be 14.5 meV 
where then energy transfer is defined as ?itj=Ei-E/. The solid 
curves are calculations using a Maxwellian distribution and a 
moderator temperature fcsr=28.8 meV. 



The factor 2 comes from the orientation factor and x" de- 
notes the susceptibility per formula unit here and in the 
body of paper. We have chosen to use the anisotropic 
form factor /^(Q) over the isotropic formiS^ because the 
anisotropic form has been carefully verified in the ordered 
antiferromagnetic system and in the Al doped YBCO sys- 
^gj^ji03ji04 Pqj, practical cross-section calculations 
is 73 mbarns sr~^. The paramagnetic description for the 
spin scattering follows from our polarized measurements 
and from our variation of the orientation factor in differ- 
ent Brillouin zones. These show no preferred polariza- 
tion of the fluctuations. The factor B{Qz) is the bilayer 
structure factor equal to sin{Qzd/2), where d is the intra- 
bilayer spacing. In the normalization to absolute units 
of the spin fluctuations below ~ 40 meV energy transfers 
we have only considered the acoustic, or odd, fluctua- 
tions to the susceptibility since these predominate at the 
momenta studied. In the notation of Fong et al. we have 
calculated x'oddi^' ^) O'^ly- 

The definitions of the magnetic cross-section used here 
are identical to those used in the phonon calibration by 
Fong et alm^. But, the definition of the q and u; inte- 
grated value of x" used by Fong et al. is larger by a factor 
of TT from the definition used in this paper and that of Dai 
et al. We have therefore removed this factor tt from the 
integral quoted by Fong et al. for the comparisons made 
hereiiS^ The vanadium calibrations conducted by Dai et 
alS^ use the isotropic form factor but otherwise the same 
equations were uscdiS^ as presented here. This allows a 
direct comparison between different experiments^Si 
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VIII. APPENDIX B: CORRECTION FOR 
HIGHER ORDER CONTAMINATION OF 
MONITOR 

The low-efficiency fission monitor in the incident beam 
only approximately removes the factor 1/fco from the 
measured cross-section. As the incident energy decreases, 
the monitor rate is perturbed because it senses an in- 
creasing fraction of higher-order neutrons reflected by 
the monochromator out of the reactor spectrum . Thus 
low-energy transfer scattering would be underestimated 
if no correction were made. To understand the relative 
intensities as a function of energy transfer, a complete 
characterization and correction for this contamination is 
essential. 

We determined the relative weight of A/2 and A/3 neu- 
trons over the energy range of our experiments from the 
relative intensities of aluminum powder Bragg peaks in 
two-axis mode. The relative flux of each component of 
the incident beam can be uniquely calculated from the 
following formulciiS^, 

I {29) cx ^(A) . . ,., E'^^(Q)'' (15) 

sm(26') sm(t^) ^ 

where $ is the incident flux, |F7v(Q)P is the nuclear 



structure factors and the sum is over constant |Q|. 

By correcting for the 1/fco monitor efficiency the rela- 
tive fraction of monitor counts in the beam from A, A/2, 
and A/3 orders were calculated as shown in Fig. We 
also computed the fractions for a Maxwellian distribu- 
tion of neutrons, a moderator temperature of kT =28.8 
meV (60° C), and for these three components only. The 
resultant fraction of each order isiSS 



2 („2)g-(n^£;.)/fcT 

F{n , E{) = ^_EjkT _|_ 4g-4£;i//cT gg-gfi./fcT ' ^^^^ 

As shown for n =1, 2, and 3 in Fig. I^Slthis simple model 
provides a reasonable description of the data given that 
we have not corrected for the energy dependent reflectiv- 
ity of the monochromator. As can be seen from Fig. 1251 
the correction factor required for low energies is substan- 
tial. Inclusion of this factor is essential for obtaining an 
accurate form for the spectral distribution, as we have 
done from to 43 meV. The correction becomes partic- 
ularly important when discussing low-energy excitations 
where it amounts to a factor two. Neglect of the factor 
would make it more difficult to observe low-energy scat- 
tering and might lead to an inference that a normal state 
spin gap was present. 
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